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ABSTRACT 

Kepler has found hundreds of Neptune-size (2-6 i?®) planet candidates within 0.5 AU of their stars. The nature 
of the vast majority of these planets is not known because their masses have not been measured. Using theoretical 
models of planet formation, evolution and structure, we explore the range of minimum plausible masses for low-density 
exo-Neptunes. We focus on highly irradiated planets with T eq > 500 K. We consider two separate formation pathways 
for low-mass planets with voluminous atmospheres of light gases: core nucleated accretion and outgassing of hydrogen 
from dissociated ices. We show that Neptune-size planets at T eq — 500 K with masses as small as a few times that of 
Earth can plausibly be formed core nucleated accretion coupled with subsequent inward migration. We also derive a 
limiting low-density mass-radius relation for rocky planets with outgassed hydrogen envelopes but no surface water. 
Rocky planets with outgassed hydrogen envelopes typically have computed radii well below 3 i?®. For both planets 
with H/He envelopes from core nucleated accretion and planets with outgassed hydrogen envelopes, we employ planet 
interior models to map the range of planet mass-envelope mass-equilibrium temperature parameter space that is 
consistent with Neptune-size planet radii. Atmospheric mass loss mediates which corners of this parameter space are 
populated by actual planets and ultimately governs the minimum plausible mass at a specified transit radius. We find 
that Kepler's 2-6 i?® planet candidates at T eq = 500-1000 K could potentially have masses < 4 M®. Although our 
quantitative results depend on several assumptions, our qualitative finding that warm Neptune-size planets can have 
masses substantially smaller than those given by interpolating the masses and radii of planets within our Solar System 
is robust. 

Subject headings: Planetary formation — Extrasolar planets — Accretion 



1. INTRODUCTION 

The first 4.5 months of Kepler data provide evidence 
for hundreds of "Neptune-size" (2-6 i?®, where i?® 
is the Eart h's radius) planets orb iting within 0.5 AU of 
their stars (Borucki et al.l l2bllal Tbh. The prevalence of 
planet candidates within this size range raises questions 
about both planetary growth and migration of Neptune- 
size planets. Assuming that many of these candidates 
are true planets, what are they, how did they form, and 
why are they so numerous? 

Kepler measures planetary sizes and orbital periods. 
In some cases, planet masses can be estimated from 
dynamical interactions between the planet and its star 
(Doppler method) or among planets in a multi-planet 
syste m (Transit Timing Variations, TTVs. Holman et al.l 
120101 ). However, the masses of most of Kepler's Neptune- 
size planet candidates will be difficult to measure. 

We model herein the growth, physical evolution, and 
interior structure of Neptune-size planets that possess 
voluminous atmospheres of light gases. Our focus is ob- 
taining estimates of the minimum plausible masses of Ke- 
pler's planet candidates. The maximum plausible mass 
of a planet of radius R p < 3i?® can be estimated from 
the mass-radius relati onship for rocky (Earth-like com- 
posit i on) planets (e.g..lValencia et aL 2006; S eager et al.l 
120071: iFortnev et ail 120071: iMarcus et al.ll2010D . In con- 
trast, estimation of minimum plausible masses requires 
more detailed modeling of planetary growth, because 
formation of low-mass planets of solar composition de- 
mands complicated and contrived scenarios involving 
large amounts of mass loss. We consider formation of 



low-density planets both through core nucleated accre- 
tion, and through outgassing of low-molecular weight at- 
mospheres. This work does not consider planet formation 
via gravitational insta bil ity ( Boss 1993) or tidal downsiz- 
ing (|NavakshiiJl20l0allbl 120111) . 

We present new core nucleated planet accreti o n cal - 
cula tions following the appro ach of IPollack et al.l ()1996f ) 
and iMovshovi tz et al.l ([2010). Whereas all previous pa- 
pers with this code emphasize the formation of Uranus 
mass and larger planets, here we present a new applica- 
tion of the code to a lower mass regime (M p < 10 M®). 
We push to lower planet masses by modeling formation 
scenarios where the gas disk dissipates well before rapid 
gas accretion. We also consider lower solid planetesi- 
mal surface densities (4 gem -2 at 5.2 AU and 6 gcm~ 2 
at 4 AU) than most previous calculations (10 gem -2 at 
5.2 AU) to attain lower heavy element core masses. Un- 
til recently, high planetesimal surface densities (about 
3 times the minimum mass solar nebula at 5.2 AU) 
were needed to model Jupiter formation on a reasonable 
time scale. Advances in the modeling of grain physics 
(jMovshovitz et aT1l2010f) now allow for a reasonable for- 
mation time for Jupiter, even with a = 4 gem -2 consid- 
ered here. 

We supplement the detailed core nucleated accre- 
tion calculations with equilibrium models of Neptune- 
size planets having H/ He envelopes ca l culated fol- 
lowing the appr o ach of iRogers fc Seagerl (|2010af ) and 
iRogers fc Seagerl ()2010bf ). The equilibrium model is less 
computationally time consuming and allows us to more 
comprehensively sample the parameter space (heavy el- 
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ement core masses, envelope masses, irradiation levels, 
and intrinsic planet luminosities) of interest. We focus 
on low-density planets having equilibrium temperatures 
of 500-1000 K, since these temperatures are relevant to 
the planet candidates found in Kepler's first quarter data 
(jBorucki et al.ll2011a| ). 

We also explore outgassing during planet formation as 
a possible origin pathway for low-density Neptune-size 
planets, and derive a limiting mass-radius (M p — R p ) re- 
lation bounding the maximum radius/minimum density 
for planets with primary de-gas sed envelopes. Following 
Elkin s-Tanton fe Seagerl J2008b), we consider outgassing 
of hydrogen gas produced when water reacts with metal- 
lic Fe in accreting materials during planet formation. 
Our outgassed exoplanet models self-consistently treat 
the connection between the planets' interior structure 
(iron core mass and mantle composition) and the mass 
of H2 degassed. We thereby provide the first exoplanet 
M p — R p relations that include the effect of an outgassed 
H2 gas layer. To derive the limiting M p — R p relation, we 
study planets that accreted from a mixture of water and 
material with chemical composition characteristic of the 
high-iron enstatite (EH) chondrite meteorite class, cor- 
responding to end member scenarios yielding maximum 
outgassed H 2 . 

We begin by describing our equilibrium model for low- 
mass planets with gas layers in Section [2] This model 
is applied in later sections to explore the mass-radius 
{Mp — R p ) relationships for low-mass planets with volu- 
minous atmospheres of light gases acquired by core nu- 
cleated accretion and outgassing of hydrogen. Section [3] 
describes the formation and properties of Neptune-size 
planets that assembled through core accretion, and Sec- 
tion |4] describes the formation and properties of planets 
that outgassed hydrogen from dissociated ices. We con- 
sider mass loss from the envelope in Section[5j We discuss 
our results and conclusions in Section [BJ 



2. MODELS OF PLANETS IN EQUILIBRIUM: 
METHODS 

We use equilibrium models — spherically symmetric 
planets in hydrostatic equilibrium — for two applications. 
The first (Section [3]) is to explore the mass-radius rela- 
tionships for low-mass planets formed via core nucleated 
accretion. The second application is to again study mass- 
radius relationships for planets that acquired an envelope 
of light gases through outgassing (Section 0} . 

Our equilibriu m model is based upon th e planet inte- 
rior model from iRogers fe S cagcr (2OlOa0). We have, 
however, included updates to the temperature profile in 
the radiative regime of the envelope, and to the outer 
boundary conditions of the planet. We use equilibrium 
models to study instantaneous states of evolving planets 
assuming the planets are undergoing quasistatic evolu- 
tion. Our work does not focus on cases where the en- 
velope dynamics or variations in the interior luminosity 
profile have an important effect. 

We assume spherically symmetric and differentiated 
planets consisting of up to four layers. From the inside 
out, these layers are an iron core, a silicate mixture, H2O, 
and a gas envelope. The coupled differential equations 
describing the mass of a spherical shell in hydrostatic 
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and the differential equation describing the radial optical 
depth, r, in the gas layer 
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are integrated inward from the top of the planet's en- 
velope. Above, m is the interior mass coordinate, r is 
the distance from the planet center, P is the pressure, 
p is the mass density, k is the mean opacity at thermal 
wavelengths, and G is the gravitational constant. 

Within each chemical layer, the equation of state 
(EOS) relates the density p{m) to the pressure P (m) 
and temperature T (m) . In analogy to the models in 
Section [21 we define the exterior boundary condition on 
the planet envelope (r = R p , m = M p , t = tr, P = Pr) 
at radial optical depth tr — 2/3. We then determine the 
corresponding pressure Pr by imposing 



Pn 



gt~R 



(4) 



where kr is the Rosseland mean opacity at the photo- 
sphere boundary, and g = GM P /R 2 is the gravitational 
acceleration. While the density p varies abruptly be- 
tween the chemical layers, both P and T are continuous 
across layer boundaries. For a specified planet compo- 
sition, energy budget, and mass, M pi the planet radius, 
R pi is iterated until a self-consistent solution satisfying 
the inner boundary condition (r — 0, m = 0) is achieved. 

We assume that the gas envelope is in radiative- 
convective equilibrium, with an outer radiative zone sur- 
rounding a convective layer at greater depths. Within 
the thin outer edge of the envelope, we adopt the 
is otropic average temperature profile from Equation (29) 
in lGuillotl (|2010f ). 
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an analytic solution to the "two-stream" gray equa- 
tions of radiative transfer for a plane-parallel irradiated 
atmosphere. The irradiation temperature, Ti rr , char- 
acterizes the short wave energy flux received by the 
planet from the star and relates through the redistri- 
bution factor, /, to the equilibrium temperature of the 
planet in the radiation field of the star T eq = / 1/,4 T,* rr 
(/ = 1/4 for full redistribution). The intrinsic tem- 
perature Ti n t = (Lp/AirRpO-B) 1 parameterizes the to- 
tal intrinsic luminosity of the planet, L p {ctb denotes 
the Stefan-Boltzmann constant). The total intrinsic 
planet luminosity, L p , is the sum total of contributions 
from envelope contraction, radioactive decay, and secu- 
lar cooling of the core. The ratio of the short-wave and 
long- wave optical depths is represented by 7. We take 
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7 = 0.6^r irr /2000 K, which |GuilloJ pOlCh found pro- 
vided a good matc h to detailed ca l culati ons of hot Jupiter 
atmospheres from lFortnev et ail (|2008l ). 

In highly irradiated planet atmospheres, the radiative 
regime of the envelope may extend to depths beyond 
where the plane parallel approximation (assumed when 
deriving Equation [5]) is valid. In these cases, once all 
of the incoming stellar radiation is absorbed at optical 
depths t 3> l/v37) we transition smoothly to the radia- 
tive diffusion equation, 
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The onset of convective instabilities 
(0 < (dp/ds) p ds/dm) determines the depth of the 
transition to the convective layer of the gas enve- 
lope. In the convective regime, we adopt an adiabatic 
temperature profile. 

We use the EOS from lSaumon et all (119951). The effect 
of uncertainties in the H/He E OS (see e.g. jMilitzer et aTl 
I2008t iNettehnann et al.ll2008l) will be small compared to 
the effect of uncertainties in the heavy element compo- 
sition and distribution, for the low-mass planets we are 
considering. While the major uncertainties in the EOS 
are at Mbar pressures or above, the pressures at the base 
of our H/He envelopes are typically less than a few tenths 
of a Mbar. As in the formation and evolution models of 
Secti on [51 we use Rosselan d mean molecular opacities 
from iFreedman et~a l. (2008). We neglect grain opacities 
in our equilibrium models, however, since we are inter- 
ested in the planet radii at late times, after all the grains 
have settled. 

Under the H/He envelope, the rock-ice interior is mod- 
eled with differentiated layers of iron, Feo.iMgo.gSiOs sil- 
icates, and H2 O. For these mater ials, we employ EOS 
data sets from Seager e t al.1 (j2007j ). which were derived 
by combining experimental data at P < 200 GPa with 
the theoretical Thomas-Fermi-Dirac (TFD) equation of 
state at high pressures, P > 10 4 GPa. The equation 
of state at intermediate pressures between ~ 200 and 
10,000 GPa is not well known, since this pressure range 
is neither easily accessible by experi ments nor by TFD 
theory. For H9 0. lSeager et all (|2007l ) used density func- 
tional theory calculations to fill in the EOS in this pres- 
sure regime, while for all other materials, they bridged 
the pressure gap by extrapolating the empirical Birch- 
Murnagham EOS and the theoretical TFD EOS to higher 
and lower press ures, respect i vely. Thermal effects are 
neglected in the iSeager et all (| 2 00 71) EOSs - at the high 
pressures found in the interior layers, thermal corrections 
have only a small effect on the density, p. An improve- 
ment over our previous models is that we now more con- 
sistently take into account the Si/Mg/Fe ratios in the 
mantle by calculating EOSs for mixtures of MgO, FeO 
and Si02- Core mass-radius relations calculated follow- 
ing this scheme (but neglecting the small contribution to 
pressure from the gaseous envelope) were also employed 
in the planet evolution calculations of Section [3] 

A major uncertainty in the validity of the models 
comes from the assumption that the layers of water and 
H (or H/He) are not mixed. For the T eq = 500 - 1000 K 
planets considered in this work, H2 and H2O are misci- 
ble at the pressures and temperatures in the model en- 



velopes. So they could, in principle, be homogeneously 
mixed. By considering the extreme where the H and 
H2O are fully separated, we set an upper bound on the 
planet radius; typically if hydrogen is mixed into the wa- 
ter la yer one expects the pl anet's radius to be smaller 
(e.g., INettehnann eTalll2010T ). Although our aim is to 
model H/He envelope planets, some of our models do 
have significant water content, and future work should 
include the miscibility of H2 and H2O. 

The planet radii in both our equilibrium and evolution 
models underestimate the planet radii measured during 
transi t in a predictable way. This "transit radius ef- 
fect" (jBaraffe et aLl2003 iBurrows et al.ll200H ) is a conse- 
quence of our exterior boundary condition (Equation |4]) , 
which pegs our model planet radii, R p at a radial optical 
depth tr = 2/3. In contrast, it is the transverse opti- 
cal depth for transmission thr ough the planet limb that 
determines the transit radius. Hansen (2008) derived a 
correction for the transit radius effect, 
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where the transit radius Rt = R p + AR is defined at a 
transverse optical depth of unity, and Hr represents the 
atmospheric scale height at the planet limb. Equation (J7J 
applies when Hr <C R P7 and assumes the outer limb of 
the planet atmosphere is well described by an ideal gas. 
For the low-mass (M p < 30 M®) planets with hydrogen- 
rich envelopes we are considering, the transit radius cor- 
rection is typically between 1% and 10%. Equation [7] 
assumes a clear cloud-free planet atmosphere; high level 
clouds or hazes could further enhance the transit radius 
effect. 

3. PLANET FORMATION BY CORE-NUCLEATED 
ACCRETION 

3.1. Methods 

Models for the formation and evolution of a planet 
consisting of a core of heavy elements and a gaseous 
envelope of solar composition are calculated acco rding 
to the procedures descri bed bv iPollack et al.l (|1996| ) and 
iMovshovitz et all (|2010f ). The formation calculation con- 
sists of three major parts: (i) the accretion rate of plan- 
etesimals onto the planet; (ii) the interaction of the in- 
falling planctesimals with the gaseous envelope; and (Hi) 
the evolution of the gaseous envelope and the determi- 
nation of the gas accretion rate. 

The planctesimal accretion rate onto the plan etary em- 
bryo is based on the equation originated by iSafronovl 
(1969). If M coro is the mass of the embryo, then the 
fundamental equation for its growth is 



dM a 



dt 



(8) 



where 7ri?^ apt is the effective geometrical capture cross- 
section, a is the surface density of solid material (plan- 
etesimals), Q is the orbital frequency, and the value of 
the g ravitational enhanceme n t fact or, F g , is obtained 
from iGreenzweig fc Lissaueri ([1992), assuming a plan- 
etesimal radius of 100 km. If no gaseous envelope is 
present, then i? ca pt = -Rcore, the heavy element core ra- 
dius. As in our previous publications, we take the feeding 
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zone from which the embryo can accrete planetesimals 
to extend 4 Hill sphere radii on either side of the orbit 
(|Karv fe Lissauerl fl994'l . and assume that the solid sur- 
face density a is constant within that zone. The value 
of a in the feeding zone is adjusted at each time step to 
take into account depletion of planetesimals by accretion 
onto the embryo and expansion of the feeding zone into 
undepleted regions, as the embryo's mass increases. 

The second element of the code calculates the inter- 
actions between planetesimals and the ga seous envelope 
of th e protoplanet as they fall through it (|Podolak et al.l 
1988). The detail s of how this c alcula ti on is performed 
are d escrib ed in iPollack et all (11996D . iHubickvi et al.l 
(|2005f) . and lMovshovitz et al.l (|2010i r These calculations 
provide the effective capture radius i? cap t to be used in 
Eq. ([8]) , as well as the deposition of mass and energy as a 
function of radius in the gaseous envelope. The effective 
capture radius can be several times larger than the actual 
solid heavy element core radius because of the effects of 
the gas on slowing down, ablating, and fragmenting the 
planetesimals. It is assumed that the material from the 
planetesimals that is deposited in the envelope later sinks 
to the heavy el ement core, releasin g gravitational energy 
in the process (IPollack et al.lll996f). This assumpt ion is 
not entirely accurate: llaroslavitz fc Podol ak (2007) show 
that some of the heavy-element material should actually 
dissolve in the envelope. Thus the 'core masses' that 
we calculate actually represent the total heavy element 
abundance in the planet in excess of the solar metal abun- 
dance; most of these heavy elements (including all of the 
rock and organic compounds) would be expected to re- 
side in the actual core of the planet. 

The third clement of the simulation is the solution 
of the four differential equations of stellar structure for 
the gaseous envelope, with energy sources from accreting 
planetesimals, from gravitational contraction, and from 
cooling. The adiabatic temperature gradient is assumed 
in convection zones. At the heavy element core bound- 
ary, the luminosity L r is set to the energy deposition rate 
for the planetesimals that hit the heavy element core. 
Outside the heavy element core, the energy supplied by 
ablated and fragmented planetesimals is included as a 
source term in the energy equation. 

At the core-envelope interface, the radius is set to that 
of the outer edge of the heavy element core. The heavy 
element core mass-radius relation is calculated using the 
equilibrium model described in Section [21 assuming 10% 
Fe, 23% Fe .iMgo. 9 Si0 3 , and 67% H 2 0, by mass. The 
heavy element core composition we adopt is motivated 
by comet compositions, and represents rock with a Fe/Si 
ratio near solar mixed with ice in a ratio of 1:2 by mass. 

At the surface, gaseous mass of solar composition is 
added at a sufficient rate to maintain an outer radius 
R p = R c ff, which is given by (|Bodenheimer et al.ll2000h 
as 
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where Rh is the Hill sphere radius, c s is the sound 
speed in the disk outside the planet, M p is the total 
planet mass, and, nominally, K = 1. Note that when 
Rh is large compared with the Bondi accretion radius, 
GMp/c^, the expression reduces to the Bondi radius, 
while in the case of the opposite limit, i? e ff — > Rh- 



In developments after the above expression was for- 
mulated, it turned out that K had to be modified. 
Three-dimensional ca lculations of disk-planet interaction 
(jLissauer et alJl2009[ ) gave the result that not all the gas 
passing through the Hill sphere is actually accreted by 
the planet; some of it simply flows through and rejoins 
the disk's azimuthal motion. The 3D simulations pro- 
vided an estimate of the effective planetary radius, which 
corresponds to K = 0.25, the value used in this paper. 

The density and temperature at the planet's surface 
are set to assumed nebular values p n cb, Tncb, respectively. 
The value of T no b is constant in time, while /9 ne b decreases 
linearly to zero with time, over a time scale Td « 2 — 3 
Myr. In a variation of this boundary condition, p no b is 
constant in time up to a time comparable to T<j, then it is 
linearly reduced to zero on a time scale of 10 5 yr. These 
assumptions roughly characterize the dissipation of the 
gaseous disk. T ne b is held constant while the planet is ac- 
creting; our model incorporates migration only through 
temperature increases subsequent to the conclusion of 
the planet's growth. Modeling simultaneous migration 
and accretion is beyond the scope of this work. 

When p no b approaches zero, the accretion of gas is 
halted and the evolution is calculated at constant mass 
over time scales up to 3-4 Gyr. The envelope mass at 
cutoff in these simulations is always small enough that 
rapid runaway gas accretion does not occur, and Eq. ((9|) 
is always valid for the determination of the gas accretion 
rate. The accretion rate required to keep R p = R g re- 
mains much lower than the limit imposed by disk physics 
in supplying materia l to the Hill sphere of the planet 
(jLissauer et al.|[2009() . Once gas accretion is shut off, R p 
rapidly falls below i? c ff , and the planet becomes isolated 
from the disk. The surface boundary condition changes 
to that of a hydrostatic atmosphere that radiates from 
the photosphere: 



L = A<kRI<t b T* s and 



kP = 



(10a) 
(10b) 



where ctb is the Stefan-Boltzmann constant, T e ff is the 
surface temperature, L is the total luminosity (energy 
radiated per second) of the planet, and k, P, and g are, 
respectively, the Rosseland mean opacity, the pressure, 
and the acceleration of gravity at the photosphere. There 
are two contributions to T e g: that from the internal lumi- 
nosity provided by the planet, and that from the energy 
absorbed from the central star and re-radiated by the 
planet. Thus: 

X, (11) 



rpi rp4 

1 off — 1 int 



where Ti nt is the internal contribution (generally small), 
and T eq is the equilibrium temperature of the planet in 
the radiation field of the star. The former quantity is 
determined from the evolutionary calculation, while the 
latter is a parameter that depends on the assumed dis- 
tance of the planet from the star and the stellar luminos- 
ity. 

The equation of state of the gas is taken from 
iSaumon et al.l (|1995l) . interpolated to our assumed com- 
position of hydrogen mass fraction X = 0.74, helium 
mass fraction Y = 0.243, and metal mass fraction Z = 
1 — X — Y = 0.017. Although the equation of state in 
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the outer, low-density layers is essentially that of an ideal 
gas, the inner regions near the heavy element core can 
be significantly non-ideal once the envelope has become 
sufficiently compressed. 

The Rosseland mean opacity calculation has three 
components. At tempe ratures above 3000 K, the mol ecu- 
lar /atomic opacities of lAlexander fc Ferguson! (|1994f ) are 
used. In practice the details of the opacities in this re- 
gion are unimportant because the energy transport is 
almost always by convection. In the temperature range 
100-3000 K the molec ular opacities, without grains, of 
iFreedman et al.l ()2008l ) are used. Grain opacities are then 
added in the temperature range 100-1800 K. Two sources 
of grains are taken into account, first, those provided by 
the ablating planetesimals as they interact with the en- 
velope, and, second, those that accrete along with the 
gas at the surface of the planet. At each time step of 
the evolutionary calculation, and at each depth in the 
envelope, the grain size distribution is recalculated, tak- 
ing into account the coagulation and settling of grains. 
The size distribution is represented by 34 bins, covering 
the size range 1.26 /im to 2.58 mm. The effective cross- 
sections for absorption and scattering are calculated as 
a function of grain size and frequency; then an integra- 
tion over grain size and frequency gives the Rosseland 
mean opacity as a functio n of depth. The details of the 
grai n physics are given in IMovshovitz fc Podo lak (2008) 
and IMovshovitz et al.l (|2010| ). The grains are composed 
purely of silicates, with a dust-to-gas ratio of about 0.01 
by mass; little error results from this assumption com- 
pared to the uncertainties in gr ain shape, sticking proba - 
bility, and radiative properties (|Movshovitz et al.1 12010). 
Grains are assumed to be completely evaporated above 
1800 K. The grains are important during the gas accre- 
tion phase. Once accretion is shut off, the grains rapidly 
settle towards the center and are evaporated. This ef- 
fect is included in the calculations and indicates that any 
grains remaining in the atmosphere have a negligible ef- 
fect upon the evolution. Thus in the final constant-mass 
evolution phase, the molecular opacities completely dom- 
inate. 

3.2. Evolution Input Parameters and Results 

The planet initially consists of a heavy element core of 
1 M® and a light element envelope of about 10~ 5 M®. 
The protoplanet is located at either 5.2 AU or 4.0 AU in a 
protoplanetary disk, with the solid surface density a = 4 
g cm~ 2 at 5.2 AU and 6 g cm" 2 at 4 AU. The initial 
evolutionary time is set to 7.3 x 10 5 yr and 4.8 x 10 5 yr, 
respectively, for a — 4 and 6 g cm -2 , approximately the 
time needed to assemble a heavy element core of mass 
M core = 1 M e . 

The quantity T nob is set to 115 K at 5.2 AU and 125 
K at 4 AU. Then p nch = axy/(2H), where a X y = 70ct 
is the surface density of the gas component. As men- 
tioned above, p no b in general declines with time. The 
scale height of the gas in the disk H = 0.05a, where a is 
the orbital distance from the star. Once started, the evo- 
lution consists of three main phases. The first involves 
primarily accretion of solids onto the heavy element core, 
with a relatively low-mass envelope and a low gas accre- 
tion rate. The solids accretion rate slows down signif- 
icantly near the point where the isolation mass (M; so ) 
for the core is reached; for a = 4 g cm~ 2 at 5.2 AU 
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Fig. 1. — Mass of the protoplanet as a function of time for Runs 
I. For Run la (black curves) the solid line denotes the mass of the 
heavy element core, the dotted line the mass of the H/He envelope, 
and the short-dash-dot line the total mass. For Run lb, the same 
line types are shown in red. 
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Fig. 2. — Mass of the protoplanet as a function of time for Runs 
II. For Run Ila (black curves) the solid line denotes the mass of the 
heavy element core, the dotted line the mass of the H/He envelope, 
and the short-dash-dot line the total mass. For Run lib, the same 
line types are shown in red. 



this mass is about 2.9 M® and for a = 6 g cm~ 2 at 4 
AU, about 2.4 M®. During the second phase, the gas 
accretion rate is about 3 times as high as the core ac- 
cretion rate, M env « 3M core , and both are nearly con- 
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TABLE 1 

Input Parameters for Evolutionary Runs 



Run 


a (AU) 


a (g/cm 2 ) 


a X Y (g/cm 2 ) 


T ncb (K) 


T d (Myr) 


M iao (M e ) 


la 


5.2 


4 


280 


115 


3.5 


2.92 


lb 


5.2 


4 


280 


115 


2.5 


2.92 


Ila 


4.0 


6 


420 


125 


2.0 


2.42 


lib 


4.0 


6 


420 


125 


0.9 


2.42 



TABLE 2 

Results from Evolutionary Runs: Masses and Radii 



Run 


M p (M ffi ) 


M corc (Me) 


M cnv (M e ) 


Rl;500 (-Re) 


-R4;500 (-Re) 


Riaooo (Re) 


#4;1000 (-Re) 


la 


8.13 


4.08 


4.05 


9.8 


8.1 


14.8 


11.6 


lb 


5.20 


3.52 


1.68 


8.0 


6.6 


15.7 


11.5 


Ila 


3.19 


2.65 


0.54 


6.0 


5.0 


17.9 


11.7 


lib 


2.66 


2.50 


0.16 


3.6 


3.3 


6.7 


6.2 



The first subscript on the radius gives the evolutionary time in Gyr. The second subscript gives the assumed equilibrium temperature of 
the planet. 



stant in time (|Pollack et al.lll~996f ). The envelope mass 
builds up relative to the heavy element core mass, which 
grows slowly. The phase of rapid gas accretion, which 
for giant planets begins once the envelope mass M env 
becomes about equal to M corc , does not occur in these 
calculations. Instead, gas accretion is cut off and the 
planet evolves through a third phase at constant mass 
with boundary conditions provided by Eq. (|10[) . During 
the early part of this phase, the planet is assumed to mi- 
grate to a position within 1 AU from the star. Represen- 
tative cases with T eq = 500 K and 1000 K are presented. 

The input parameters of the four runs are shown in 
Table 1, which includes er, the gas surface density uxy, 
the surface boundary temperature T ne b, and the isolation 
mass M- lso . 

The results of our calculations for Runs la, lb, Ila, and 
lib are shown in Figures [TJ [2j and [3] The masses and 
radii that are derived for the four runs are listed in Table 
2. 

Run la is based on a disk with a = 4 g cm -2 at 5.2 AU. 
This value is only slightly greater than that of the min- 
imum mass solar nebula. But note that our calculation 
of F g (Eq. [8|) neglects transport of solids into or away 
from the planet's accretion zone. Moreover, our planetes- 
imals are all assumed to have the same radius, 100 km. 
In fact there must be a range of planetesimal sizes, and 
the effective planetesimal size is not well known. Smaller 
planetesimals would result in mor e rapid accretion (see 
footnote 3 of iLissauer et al.l I2009T ) . The accretion rate 
that is actually calculated may thus correspond to a value 
of a slightly different from 4 g cm~ 2 . 

The details of the calc ulation with the parame ters of 
Run la are presented in iMovshovitz et alT (|2010f ) , their 
Run cr4. In that paper the run is continued well into 
the phase of rapid gas accretion, and is terminated with 
M core = 4.74 M® and M env = 34 M®. The formation 
time for a giant planet is found to be 4 Myr. In the 
present run, the accretion of gas and solids is cut off at 
a time of 3.5 Myr, co nsistent with estim ated lifetimes of 
protoplanetary disks (jHillenbrandl 120081 ) . At that time 
the value of p ne b is assumed to decrease to zero on a 
time scale of 10 5 yr. The calculation is then continued up 
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Fig. 3. — The protoplanet's total luminosity, including internal 
and irradiation contributions, as a function of time during the for- 
mation phase and the contraction/cooling phase for Run la (solid 
curve), and Run Ila (dashed curve). The equilibrium temperature 
is increased to 500 K, after the formation phase, during these runs. 



to Gyr times with constant values of M corc = 4.08 M® 
and M onv = 4.05 M®. At the beginning of this phase 
the equilibrium temperature is gradually increased, on a 
time scale of 4 Myr, to an assumed final value of 500 K. 
A gradual increase in T oq to 1000 K was accomplished in 
a total time of 6 x 10 7 yr. The final values of R p for these 
two temperatures and for times of 1 and 4 Gyr are given 
in Table 2; they are close to Jupiter's radius R,/ w 11 i?®, 
even though the planet's mass is only 8.13 M®. 

Run lb also is based on the run cr4 from 

IMovshovitz et all (|2010D . In this case the accretion of 
gas and solids was cut off at 2.5 Myr, at which point 
M core = 3.52 M® and M env = 1.68 M®. The evolution 
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was again continued into the phase of cooling and con- 
traction at constant mass, with assumed values of T eq of 
500 and 1000 K. In the case with T cq = 1000 K, the final 
radii are again comparable to or larger than Rj. In the 
case with T oq = 500 K the minimum radius is 6.6 i£®, 
only slightly smaller than the corresponding value in Run 
la. 

Run Ha is an entirely new calculation, with the planet 
forming at 4 AU in a disk with a — 6 g cm~ 2 . During 
the initial phase of rapid core accretion, the luminosity 
reaches a maximum of 3.1 x 10~ 7 L© at a time of 6.2 x 10 5 
yr. The heavy element core mass is 2.2 M© and the 
core accretion rate M cor c = 5 x 10~ 6 M© yr -1 at this 
time. Later, at 1 Myr, M core has decreased to 2 x 10~ 7 
M© yr -1 , and M env has increased to 5 x 10~ 7 M© yr -1 . 
The luminosity has decreased to 10~ 8 L©. Because of 
computational time limitations, and to obtain a lower 
envelope mass than that found for Run lb, the accretion 
in this run was cut off at 2 Myr, with M core = 2.65 
M© and M env = 0.54 M®. If the evolution had been 
continued up to 2.5 Myr, the heavy element core mass 
would have been practically unchanged, and the M env 
would have increased by about 0.25 M©. At the end of 
the contraction/cooling phase the radii are in the range 
5-6 i? e for the case of T eq = 500 K, and for T cq = 1000 K 
they are larger then Rj, close to the values obtained in 
Runs I for that temperature. 

To investigate the effect of an even smaller value of 
M cnv , Run lib was calculated with the same parameters 
as Run Ha, but with an arbitrary accretion cutoff at 9.1 x 
10 5 yr. At that point, M core = 2.5 M© and M cnv = 0.16 
M®. Final radii turned out to be in the range 3-4 i?© 
for T eq = 500 K and in the range 6-7 R B for T oq = 1000 
K. The significant reduction in envelope mass resulted 
in final radii that are about half the values obtained for 
Run Ha. 

We neglect heating from radioactive decay in the core 
nucleated accretion calculations. Including this addi- 
tional energy source would delay envelope contraction 
and planet cooling. Consequently, the planet radii at 
1 Gyr and 4 Gyr in Table 2 may be systematically un- 
derestimated by a small amount. We estimate that, for 
the cases in Table 2, the planet luminosity from radioac- 
tive decay would be roughly one order of magnitude 
smaller than luminosity from envelope contraction, as- 
suming bulk Earth a bundances of K. U , and Th in the 
heavy element cores (|Van Schmuslll995l ). The fractional 
contribution to the planet energy budget from radioac- 
tive heating will be higher for older planets (4 Gyr) and 
cases where the heavy element core contributes a larger 
fraction of the planet mass (Run II). 

3.3. Equilibrium Model Results 

In this section we explore planet radii over a wide 
range of heavy element core masses, envelope masses, 
irradiation levels, and intrinsic planet luminosities. The 
planet formation and evolution model described in the 
Section 13.11 is computationally intensive. Since it is not 
feasible to simulate planets under all conditions of inter- 
est following that approach alone, we enlist an equilib- 
rium planet structure model (Section [5]) to cover a wider 
range of parameter space. 

Our equilibrium model shows good agreement with the 



planet evolution models in Section 15721 despite the differ- 
ences in their treatment of the outer radiative regime, 
the intrinsic planet luminosity, and the effects of stellar 
insolation. For each entry in Table [5J we applied the 
equilibrium model to simulate the same combination of 
M core , M env , T eq and T int . The radii at T eq = 500 K in 
the two models agree to better than 0.2 i?© in every case. 
The planet radii at T eq = 1000 K are more sensitive to 
model assumptions and exhibit larger discrepancies (up 
to 14%, with the equilibrium model radii systematically 
below those in Table [5]) . 

We explored the parameter space of M core , M env , T eqjL 
and Ti n t with our equilibrium model. Figures U and |o] 
present a selection of mass-radius (M p — R p ) curves at 
(a) T eq = 500 K and (b) T eq = 1000 K. Figure H displays 
the effect on the radius of varying the envelope mass 
fraction, while Figure [5] shows the effect of varying the 
planet's intrinsic luminosity, L p ~ AirR^aT^. The thick 
solid line is common between FigureU]and[Sl representing 
M env = 0.2M p and L p /M p = lO" 105 Wkg" 1 . L p /M p = 
10~ 10 ' 5 Wkg -1 corresponds to both the 8.3 M© evolu- 
tion model (Run la) at 4 Gyr, and the 3.19 M© evolution 
model (Run Ha) at 1 Gyr (independent of T eq ). 

The M p — R p curves for low-mass planets with volumi- 
nous gas layers show several notable features. First, the 
planet radii (at constant envelope mass fraction, T eq , and 
L p /M p ) increase dramatically toward low planet masses. 
This is due to the low surface gravities, and thus large 
atmospheric scale heights found at low masses. Second, 
the radius of planets having identical envelope mass frac- 
tions, M env /M p , are remarkably insensitive to the planet 
mass when M p > 15 M©. At these masses, increased 
compression of the envelope offsets the effect of increas- 
ing the planet mass. Third, for planets of identical total 
mass (within the mass range plotted) the planet radius 
increases monotonically with the envelope mass fraction. 
Fourth, Ti n t and T eq both have a stronger effect on the 
radius of low mass planets compared to their more mas- 
sive counterparts. This is understandable, because, given 
the same envelope mass fraction, in lower mass planets 
the envelope accounts for a larger fraction of the planet 
radius. 

Planet radii between 2 and 6 i?© are of special interest, 
because Kepler is finding a l arge number of planet candi- 
dates within this size-range (|Borucki et al.ll2011al lbT) . We 
plot in Figure [6] combinations of M env and M p that yield 
planet radii within this range. Planets at 2 i?© can con- 
tain at most 0.08% of their mass in H/He at T eq = 500 K, 
and at most 0.0015% at T eq = 1000 K. Larger planets 
can support more massive envelopes. A 6 i?© planet at 
T eq — 500 K requires an envelope accounting for at least 
a few percent of the planet mass. At T eq = 1000 K and 
6 i?©, between 0.1% and 23% H/He by mass is possible, 
depending on the planet mass and intrinsic luminosity. 

It is important to note that the M p — R p relations in 
Figures 31 and [5] are not isochrons, but correspond 
instead to constant total intrinsic luminosity per unit 
mass, Lp/Mp. The total intrinsic luminosity, L p , is the 
sum total of heating from radioactive decay, cooling of 
the planet core, and contraction of the planet envelope. 
In the evolution calculations from Table [21 the planet 
luminosity contribution from envelope contraction alone 
ranges from 10 -9 8 to 10~ 112 Wkg -1 at 1 Gyr and from 
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Fig. 4. — Equilibrium mass-radius relations for various choices of envelope mass fraction, M en v /Mp. All data in this plot have L p /M p = 
10 -10.5 Wkg-\ and (a) T eq = 500 K or (b) T eq = 1000 K. Each curve corresponds to a different value of M env /M p : 0.001 (thin solid), 0.01 
(thin dashed), 0.05 (thin dot-dashed), 0.1 (thin dotted), 0.2 (thick solid), 0.3 (thick dashed), 0.4 (thick dot-dashed), and 0.5 (thick dotted). 
Black lines denote our model radii (defined at a radial optical depth r = 2/3), while the corresponding blue lines represent radii corrected 
for the transit radius effect. The thick red line is the mass-radius relation for icy heavy element cores having no envelope (M e nv = 0). Red 
triangles present the subset of Table [2] evolutionary run results that have L p /M p rs 10 -10 5 Wkg -1 : Run la (M p = 8.3 A%) at 4 Gyr, 
and Run Ila (M p = 3.19 Mgj) at 1 Gyr. The green curves show the effective planet Roche-lobe radius for four different choices of host-star 
properties representative of spectral classes M5 V, M0 5V, K0 V, and G2 V (in order of increasing Roche lobe radii). The K0 V and G2 V 
Roche-lobe radii are beyond the scale of the T eq = 500 K plot. 
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Fig. 5. — Equilibrium mass-radius relations for various choices of intrinsic planet luminosity L p /M p 
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1(T 10 - 5 to 1(T 12 - 4 Wkg" 1 at 4 Gyr. Some of the low- 
L p curves in Figures [5] and [6] do not extend to higher 
masses because they encounter unphysically low planet 
interior entropies. Although L p is a proxy for the age of 
the planet, the relationship between L p and planet age 
depends on the planet's mass, composition, abundance 
of radioactive isotopes, insolation history and dynamical 
history. Since our equilibrium models are presented at a 
specified L p , we have side-stepped the issue of relating 
L p to planet age and present the model radii in a way 
such that they can be applied to many different evolution 
scenarios. Our aim with the equilibrium models is to 



broadly explore parameter space; it is beyond the scope 
of this work to relate L p and age directly by simulating 
all possible planet evolution histories. 

Simulated planet radii for planets at T eq = 1000 K 
may be in error by up to 20%. The problem is in extrap- 
olating the opacity tables at the high pressure end. This 
in turn makes the radiative-convective boundary uncer- 
tain (a deeper radiative-convective boundary makes for a 
smaller planet). Planets at T eq — 500 K are less affected 
by this opacity-caused radius problem (< 10% radius un- 
certainty for M p > 3 M©). This issue affects both our 
equilibrium and evolution models. 
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Fig. 6. — Planet mass and envelope mass that are consistent with 
a particular planet radius, for planets comprised of ice-rock inte- 
riors surrounded by H2 and He in protosolar proportions. These 
models represent planets that formed beyond the snow line by core 
nucleated accretion. We plot the envelope mass fraction as a func- 
tion of total planet mass for planets with radii (a) R p = 2 , (b) 
4 , and (c) 6 -R® . Black curves represent planets at T eq = 500 K, 
while red curves correspond to T eq = 1000 K. The line style indi- 
cates the planet luminosity: L p /M p = 10 -11 Wkg -1 (dashed), 
L p /Mp = 10" 10 Wkg" 1 (dot-dash), and L p /M p = 10~ 9 Wkg- 1 
(solid). The thin dotted lines are contours of constant envelope 
mass loss timescale, = M env /M. Each contour is labeled with 
log(t Af /Gyr) for tLxuv I Lbol = 10~ 6 , and can easily be scaled 
for other choices of tLxuv / Lbol using Equation (1 1 3 D ■ 



PLANET FORMATION BY OUTGASSING OF 
HYDROGEN 

4.1. Model 



Outgassing provides a mechanism for low-mass terres- 
trial planets to acquire an atmosphere even if they fail to 
accrete H and He from the protoplanetary nebula. In this 
section we explore the optimum conditions for a planet 
to acquire a voluminous gas enve lope through outgassing. 
We bas e our model approach on Elkins-Tanton fe Seagerl 
(2008b. a), with the improvements of a more detailed in- 
terior structure model and a calculation of the planet 
radius. 

We focus on outgassing of H2 produced when wa- 
ter reacts with meta llic Fe in accreting materials dur- 
ing planet formation (iRi ngwood 19791: fWanke fe Dreibusl 
U994tlElkins-Tanton fe Seagerll2008bl ). Hydrogen gas has 
the potential to yield the most voluminous outgassed at- 
mospheres, being both of low-molecular weight and (for 
some planetesimal compositions) degassed in substantial 
quantities. Although we do not consider these processes 
in detail here, in general, outgassing may also proceed 
during accretion as impinging planetesimals are heated 
and vaporized upon impact; during magma ocean solid- 
ification as volatiles are partitioned between the atmo- 
sphere and melt; and during volcanic/tectonic activity 
after the planet has formed. 

The reaction between water and metals during plan- 
etary accretion and differentiation intrinsically links the 
planet's interior structure to its initial atmosphere's mass 
and composition. Metallic iron forming the planet will 
either differentiate to contribute to the planet iron core, 
or become oxidized and incorporate into the planet man- 
tle. Given an initial composition for the primordial ma- 
terial forming a planet, there are two extremes to the 
eventual planet outcomes. If none of the available wa- 
ter and metals in the accreting materials react (reduc- 
ing conditions), the planet will have a maximally mas- 
sive metallic core, relatively iron-poor mantle, minimal 
outgassed H2, and maximal leftover H2O. In contrast, if 
the water and metals react to the maximal extent pos- 
sible (oxidizing conditions), the planet will have a mini- 
mal iron core mass, iron-rich mantle, maximal outgassed 
H2, and minimal leftover H2O. When Fe is the limiting 
reagent, this extreme will corres pond to a coreless planet 
(|Elkins-Tanton fe Seagerll2008aft . 

To bound the radii of outgassed rocky planets, we con- 
sider the end-member case of a planet formed purely 
from high iron enstatite (EH) primordial material. The 
motivation for this choice is three- fold. First, out of 
all meteoritic compositions, EH material has the poten- 
tial to degas the most H 2 per unit mass (up to 3.6%, 
lElkins-Tanton fe Seagerll2008bl) Second, the oxygen iso- 
tope mixing model (jLodders! [2000) predicts that the 
Earth accreted from material that was 70% EH chon- 
dritic matter by mass. Third, heating of EH mate- 
rial releases a low mean molecular weight atmosphere; 
ISchaefer fc l£glevl ([20101 ) calculated 44% H 2 , 31% CO, 
17% H 2 0, 5% CO2, and 3% other molecules by volume 
under their nominal conditions (1500 K, 100 bars). Thus, 
complete oxidation of an EH planet should achieve ef- 
fectively the maximum radius plausible for planets with 
outgassed atmospheres. 

For the EH material we adopt the chemical c o mpos i- 
tion of meteorite ALHA77295 from Uarosewichl (| 19901) . 
We distill the mineralogy in our model to include only 
the most plentiful and important constituents: metallic 
Fe, FeS, FeO, Fe 2 3 , MgO, Si0 2 , H 2 0, and H 2 . Fol- 
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lowing an approach similar to Sot in et all (j2007|), less 
abundant elements are represented by their most similar 
neighbors in the periodic table: metallic Ni is added to 
metallic Fe, Ca is added to Mg, and Al is divided equally 
(by number) between Si and Mg to preserve charge con- 
servation. Other trace constituents (TiC>2, Cr2C>3, MnO, 
Na20, K2O, P2O5, Co, which combined account for less 
than 2.2% by mass) are neglected. The resulting sim- 
plified composition adopted for the primordial rocky EH 
planctesimals consists of (by mass) 38.2% Si02, 25.2% 
metallic Fe, 14.3% FeS, 20.6% MgO, 1.7% H 2 0. Note 
that H2O included in the EH material is adsorbed to the 
surface or chemically bound to the minerals. 

We consider planets initially formed from a mixture 
of EH material and H2O ice. The H2O ice is in ad- 
dition to the 1.7% H2O by mass included in the EH 
minerals. We compute the planet bulk composition af- 
ter outgassing from stoichiometry (Table |3|), assuming 
some fraction of the accreted iron reacted with water 
(Fe + H2O — >■ FeO + H2) before sinking to form the 
planet's metallic core. We note that although we con- 
sider only Fe in our reduced EH chemical composition, 
Ni can also form oxides and be incorporated in silicates. 
Nickel accounts for 8% of the generalized metallic Fe in 
our distilled EH chemical composition - the Ni abun- 
dance is abundance in ALHA77295 is 1.83% by mass. We 
do not vary the S mass fraction of the iron core in our 
models, effectively assuming metallic Fe and FeS oxidize 
in equal proportions. We do not follow any S released in 
the conversion of FeS to FeO. 

Our interior models of outgassed planets comprise up 
to four chemically distinct layers: an Fe/FeS core, silicate 
mantle, water layer, and hydrogen atmosphere. The bulk 
chemical composition of the planet after outgassing de- 
termines the relative masses of the planet layers and the 
composition of the silicate mantle. All of the degassed 
H2 is included in a gas layer surrounding the planet. We 
place all of the FeS and metallic iron in the planet core. 
We model the H2O in a differentiated water layer sur- 
rounding the mantle, although in practi ce some water 
may be sequestered into the silicates fe.g.. lElkins-Tantonl 
200-x and references therein). All of the remaining 
species (Si02, MgO, FeO, Fe20s) make up the mantle. 
The ratio of MgO /FeO sets the Mg # of the silicates (Mg 
# = Mg/(Mg+Fe) by number). We adjust the mantle 
equation of state to reflect the relative abundances of 
Si02, MgO, FeO and Fe2 03, modeling the silicates as 
a mix ture of (Mg,Fe)0 magnesiowustite (lElkins-Tantonl 
I2008D . Fe 2 3 h ematite (IWilburn fc Bassettl I1978D. and 
stishovite Si0 2 (jAndrault et al.l 11998ft . Outgassed bulk 
compositions and the corresponding planet properties are 
reported in Tabled 

4.2. Results 

We find that planets accreted from solid bodies that 
were abundant in our solar nebula can degas at most 
1.7% of their mass in H2. This limit obtains for a fully- 
degassed coreless EH composition planet that accreted 
just enough additional water (13.0% by mass) to fully 
oxidize all available iron to Fe203. EH material alone 
does not contain sufficient H2O on its own to oxidize all 
the metallic Fe within its bulk (only up to 15.2% of the 
Fe). The accreted material must include an additional 
8.6% H2O by mass in order to convert all the metallic 



Fe into FeO, or an additional 13.0% H 2 by mass to 
convert all the metallic Fe into Fe203. With any more 
water than this, the metallic Fe becomes the limiting 
reagent. The maximal outgassed H2 atmosphere that we 
derive he re is slightly lower than the value 3.6 wt % H2 
found bv lElkins-Tanton fe Seagerl (|2008bf ). Differences 
in the representative EH chemical compositions assumed 
account for this disparity. 

Mass-radius relations for planets harboring H2- 
envelopes from outgassing are shown in Figure [7] at both 
T eq = 500 K (Figure EK) and T eq = 1000 K (Figure [7b). 
The blue dot-dashed curve provides an upper limit on the 
radius of planets accreted from primordial chondritic ma- 
terial alone (without additional water ice), corresponding 
to the extreme where the oxidizing reaction proceeds un- 
til all of the H2O bound to the minerals is expended and 
0.2% of the planet mass is released in H2. After accret- 
ing enough additional water (13% by mass) to convert 
all available Fe to Fe203, the magenta solid line repre- 
sents planets having the maximal fraction of their mass 
(1.7%) in a degassed H2 envelope. This curve may be 
taken to bound the maximum radius/minimum density 
relation for planets with de-gassed H2 envelopes, but no 
freeH 2 0. 

Planets that accreted more than 13.0% by mass wa- 
ter with the EH chondrite material would have water 
left over even if all the metals in the planet iron core 
were expended in the outgassing reaction. In Figure [7] 
we show M p — R p relations of an example with initially 
20% by mass water ice in the primordial composition 
(dotted curves). The fully degassed planets with excess 
water have, in fact, a lower average density compared to 
the planets with the highest mass fraction of de-gassed 
H2 - the effect of the lower density ice-rock interior off- 
sets the decreased proportion of H2. In Figure we 
model the H2O layer as a distinct chemical layer below 
the outgassed H 2 envelope, but mixing of H 2 and H 2 is 
another possibility. If H2O and H2 are mixed in the en- 
velope, the planet radii would be smaller than the model 
radii in Figure [7] due to the decreased atmospheric scale 
height compared to the differentiated case. 

The radii of the outgassed planets depend on the 
intrinsic luminosity of the planet. In Figure we 
show mass-radius relations for planets with L p /M p = 
^Q-10.5 Wkg -1 . Increasing (decreasing) the planet's in- 
trinsic luminosity by a factor of 10 affects the planet 
radii in Figure[7jby at most +16% (-9.5%) at 5 M© and 
+4.5% (-3.2%) at 30 M®. Small planet masses and high 
H2 contents both increase the radius dependence on L p . 

We explore in Figure [5] the mass of H2 required by EH 
composition planets to reach radii of 2 to 3 i?®. Figure [8] 
is the outgassing analog to Figure[6]for core nucleated ac- 
cretion. In Figure we restrict our attention to planets 
without significant amounts of H2O on their surface or in 
their envelopes. The envelope mass fractions, M env /M p , 
at a specified radius are not strongly sensitive to the dis- 
tribution of Fe within the planet interior (i.e., whether 
the Fe is differentiated in the metallic core, or included 
in the mantle as oxides) - we show the case where all Fe 
is oxidized to FeO. Upper bounds on the H 2 wt % for 
several of the limiting cases in Table [3] are indicated by 
colored horizontal lines. 

Our main conclusion from this section is that planets 
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TABLE 3 

Bulk compositions of EH-composition planets with outgassed H2 envelopes. The first column represents the fraction of 

ACCRETED IRON THAT IS OXIDIZED AND INCORPORATED IN THE PLANET'S MANTLE. THE NEXT FOUR COLUMNS GIVE THE COMPOSITION OF 
THE EH PLANET AFTER OUTGASSING, ASSUMING ALL OUTGASSED H2 IS RETAINED. NEGATIVE ENTRIES IN THE H2O COLUMN INDICATE A 
WATER DEFICIT AND REPRESENT THE PROPORTION OF ADDITIONAL WATER (BEYOND WHAT IS INCLUDED IN THE EH MATERIAL) THAT NEEDS 
TO BE ACCRETED IN ORDER TO OXIDIZE THE SPECIFIED FRACTION OF IRON. THE LAST FIVE COLUMNS REPRESENT THE CHEMICAL MAKE-UP 
OF THE SILICATE MANTLE, AND DETERMINE THE MANTLE EQUATION OF STATE. In ROWS 2-4 WE NEGLECT FE2O3 AND ASSUME ONLY FeO IS 
PRODUCED WHEN IRON IS OXIDIZED. ROW 5 LISTS THE EXTREME END-MEMBER CASE WHERE ALL IRON IS OXIDIZED TO FE2O3. 
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Fig. 7. — Mass-radius relations for exoplanets with outgassed H2 envelopes. The planets are assumed to have formed purely from a 
combination of EH chondrite material and water ice. Accreting material with 20% water ice by mass (dotted lines), 13% water ice by mass 
(solid lines) 8.6% water ice by mass (dashed lines), and no additional water ice (dot-dashed lines) are considered. The line color indicates 
the fraction of accreted iron that reacted with water. Black corresponds to planets with no outgassed H2 and a maximally massive iron 
core (0% Fe reacted). Blue corresponds to planets where 15.2% of the Fe reacted - the maximum amount possible for pure EH material 
without added water. Green represents an end-member case wherein all the metallic Fe that accreted to the planet is converted to FeO. 
Finally, magenta lines correspond to planets that outgassed the maximum possible H2 for their initial chemical makeup - 100% of their 
accreted iron is oxidized to Fe2 03. Both the green and magenta M p — R p relations represent core- less planets, but they differ in the 
oxidation state of iron inside the planet (FeO versus Fe20a) and in the overall proportion of H2 released. Planet equilibrium temperatures 
of a) T eq = 500 K and b) T eq = 1000 K are shown. A fiducial intrinsic luminosity L p /M p = 10~ 10 ' 5 Wkg -1 is assumed in all cases. These 
curves do not include atmospheric escape of H2 . 



of mass M p < 30 M e with outgassed H atmospheres 
typically have radii less than 3 i?© (Figures [7] and [8] ) . 
Larger radii are found at the low-mass extreme of the 
M p — R p relations in Figure [71 but correspond to planets 
with very tenuous, loosely bound, envelopes. Outgassing 
of H2 from planets accreted from rocky material alone 
most likely cannot account for the Kepler planet candi- 
dates with radii between 3 and 6 i?®. 

5. MASS LOSS FROM LOW-DENSITY ENVELOPES 

A major question is whether the high T eq , light ele- 
ment, low gravitational binding energy envelopes mod- 
eled above are stable and could be retained over gigayear 
timescales. It is precisely in the low-mass, low-molecular 
weight, high T eq regime we are considering in which plan- 
ets are expected to be most susceptible to mass loss. Be- 
low we consider, in turn, the importance of Roche-lobe 
overflow, and XUV-driven atmospheric escape. 

Roche-lobe overflow can limit the radii of low-density 



planets at close orbital separations from their host stars. 
Our planet interior model assumes spherical symmetry 
and neglects tidal forces, but this approximation starts 
to break down for planets near their star. The effective 
radius of a planet's Roche lobe is approximated by 



0.49g 



2/3 



a 0.6g 2 / 3 + In (l + q 1 / 3 ) 



0.49g 1/3 - 0.04V /3 



(12) 

where q = M p /M* (Eg gletonllT98l . The Roche-lobe ra- 
dius sets a firm upper limit on the planet radius; any 
material outside the planet's Roche lobe is not gravi- 
tationally bound to the planet and can escape. We plot 
planet Roche-lobe radii in Figures|4]and [5]for a sampling 
of representative host star properties: G2 (1 M Q , 1 £©), 
K0 (0.79 M , 0.552 L L M0 (0.51 Mm 0.077 £©), M5 
(0.21 M , 0.0076 L e ) (ICarroll fe Ostliei 120071) . In com- 
puting the Roche-lobe radii, we have assumed a plane- 
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Fig. 8. — Planet mass and outgassed H2 envelope mass that 
are consistent with a particular planet radius, for EH-composition 
planets without H2O on their surface or in their envelopes. We plot 
the envelope mass fraction as a function of the total mass of the 
planet for planets with radii (a) R p = 2 Rq, (b) 2.5 R®, and (c) 
3 Rq. Horizontal lines indicate the maximal H2 wt % degassed in 
three limiting cases: if all H2O adsorbed in the EH material reacts 
with metals (0.2%, blue), if all Fe in the EH material is converted 
to FeO (1.2%, green), and if all Fe in the EH material is converted 
to Fe203 (1.7%, magenta). This figure is the outgassing analog 
to Figure ED for core nucleated accretion, and all the red and black 
lines follow the same naming conventions. Black curves represent 



planets at T e , 



500 K, while red curves correspond to T e , 



1000 K. The line style indicates the planet luminosity: L p /M p = 
10" 11 Wkg" 1 (dashed), L p /M p = 10~ 10 \Vkg~ 1 (dot-dash), and 
L p /M p = 10 — 9 Wkg -1 (solid). The thin dotted lines are contours 



of constant envelope mass loss timescale, t. 



M env /M. 



contour is labeled with log(f.,v/Gyr) given eLxuv I Lbol 



Each 
10~ 6 . 



tary albedo A — when relating T eq to the semi-major 
axis, a; taking reflection into account with A ^ will re- 
sult in smaller semi-major axes and smaller tl- Roche- 
lobe overflow is not an issue for T eq = 500 — 1000 K 
planets surrounding a solar analog star. In contrast, 
when orbiting an M star many of our low-density low- 
mass planets do fill their Roche lobes. Our equilibrium 
planet models are not a priori pegged to a given star 
spectral type. Tidal effects and the Roche-lobe radius 
set a lower bound on M* for which the low M p tail of 
our equilibrium models are applicable. 

XUV-driven mass loss is expected to be very impor- 
tant for low mass, low density planets. This results from 
the combined effect of large cross-sections to stellar ir- 
radiation, low surface gravities, and low envelope bind- 
ing energies. Predictions for the exoplanet mass loss 
rates suffer from unknowns in the stellar XUV fluxes, 
the conditions at the planet exosphere, and the mass 
loss e fficiency. We consider energy-limited mass loss , 
(e.g.. iLammer et all 120031: iLecavelier Pes Etangsl l2007t 
Vale ncia et al.ll2010D , 



M = - 



eirFxuvRxuvRp 



GMpKt 



pJ^-tide 



(13) 



The efficiency e represents the fraction of the energy 
in XUV photons incident on the planet that goes into 
unbinding particles in the planet atmosphere; we take 
e = 0.1, but M can easily be rescaled to another choice 
of e. Fxuv represents the flux of photoionizing radiation 
impinging on the planet. K tl de is a correction factor that 
accounts for tidal effects in the Roche potential of plan- 
ets in close proximity to their star (given by equation 
17 in lErkaev et aLI 120071 ). Finally, R X uv reflects the 
planet radius at which XUV photons are absorbed. We 
estimate Rxuv following order-of-magnitude argum ents 
gleaned from Section 2 of iMurrav-Clav et al.l (|2009D . 
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uv 
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PrRxuv 
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(14) 



where Njj ~ 5 x 10 21 mT 2 is roughly the column of 
neutral hydrogen needed to reach rxuv ~ 1, Pr is the 
pressure at R p , and Hr is the pressure scale height at 
R p (where r ~ 1 for visible light). 

We take an illustrative example of planets orbiting a 
solar analog star to explore the order of magnitude of 
mass loss rates. Figure [5] shows estimated mass loss 
rates for the planet models presented in Figure |4] For 
our assumed solar- twin host star, we compute Fxuv for 
T eq by scaling the integrated solar XUV flux measured 
bv lRibas et all ((2005T) (F X uve = 4 -6 x 10~ 3 Win" 2 at 
1 AU). We find that, for L XU v/Lbol = 3.4 x 10~ 6 = 
Lxuvq/Lbolq and e = 0.1, planets at the low mass 
extreme of our M p — R p relations have implausibly short 
envelope mass-loss timescales t,v = M env /M < 1 Gyr. 

We use energy limited mass loss (Equation [T3")) to in- 
clude contours of constant log(i Af /Gyr) in Figures [6] 
and[U The contour values represent log (t^/ Gyr) cor- 
responding to (eLxuv I Lbol = 10 -6 ), but can easily be 
scaled to reflect other parameter choices: 



tjfo oc {eLxuv I L 



BOL 
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At a specified T eq , the t^j contours are independent of 
the host star mass so long as tidal effects can be neglected 
{Kude ~ 1)- For the (T eq , R p ) combinations sampled in 
Figure |6] and this approximation holds for main se- 
quence host stars that are KO V or earlier, but breaks 
down for M stars. We emphasize that gives an in- 
stantaneous measure of the time that the planet would 
take to lose its envelope at the calculated current mass 
loss rate. M is expected to vary over a planet's lifetime. 
Stars that are more active (e.g., younger) than our Sun 
would have higher photoionizing fluxes. 

We find that planets at the low mass extremes of Fig- 
ures |U and |U have short envelope mass-loss timescales 
= M env /M < 1 Gyr (assuming eL X uv I 'Lbol = 
10~ 6 ). One could conceivably choose a threshold enve- 
lope loss timescale and then derive a lower bound on 
the planet mass at a given radius based on that assump- 
tion. We elaborate this possibility further in Section [Ol 

6. DISCUSSION 
6.1. Formation of Low- Density Neptune-Size Planets 
6.1.1. Core Nucleated Accretion 

Can core nucleated accretion form low-density plan- 
ets in the size range of 2-6 -R©? The answer is yes, 
given appropriate conditions. The solids surface den- 
sity in the protoplanetary disk must be appropriate for 
the accretion of heavy element cores a few times as mas- 
sive as Earth. These cores must grow early enough to 
accrete significant gaseous envelopes, but gas accretion 
must end early enough to avoid runaway gas accretion. 
Our evolution calculations in Section [3] demonstrate that 
M p < 10 M© H/Hc rich planets can form for plausible 
choices of a and disk lifetimes. The values we chose for 
a are only slightly above that in the minimum-mass so- 
lar nebula, but high enough so that Jupiter at 5 AU can 
form in 4 Myr. 

A second related question is whether core nucleated ac- 
cretion with subsequent migration can lead to Neptune- 
size planets at high irradiation temperatures T eq > 
500 K. Our evolution calculations uncover two factors 
that complicate achieving 2-6 i?© planets following in- 
ward migration. First, high irradiation temperatures 
lead to very large fluffy planets with radii R p > 6 i?©. 
Second, very long migration timescales are required to 
heat a planet to T eq = 1000 K while keeping its envelope 
intact. We elaborate both of these points below. 

The salient feature of our evolution calculations is that, 
despite M p < 10 M®, the irradiated planet radii at 1 and 
4 Gyr are, in many cases, larger than 6 R<§. Specifically, 
all cases in Table H with T eq = 1000 K or M p > 5 M© 
have radii in excess of 6 i?©. Lower mass envelopes are 
required to yield Neptune-size planets at these high irra- 
diation levels (Figure [6]) . Truncating gas accretion ear- 
lier (shorter disk lifetime) and subsequent envelope mass 
loss are two potential avenues toward R p < 6 i?© . While 
the model radii at 1000 K are very uncertain due to un- 
certainties in the opacities near the radiative-convective 
boundary, for the cases in Table [5] the conclusion that 
R p > 6 i?© is, nonetheless, robust. 

We found that slow planetary migration is needed for 



the low-mass envelopes to stay bound as the temperature 
at the planetary surface increases. In our evolution calcu- 
lations, the planets initially assemble at T ne \> = 115 K or 
125 K and then migrate inward to T eq = 500 K or 1000 K. 
The long migration timescale (~ 40 Myr) taken to reach 
T eq — 1000 K with the envelope intact is in tension 
with typical disk lifetimes (1-10 Myr). The migration 
timescale to reach 500 K (~ 5 Myr) is more plausible. 
It is possible that evaporative cooling or increases in the 
envelope mean molecular weight from preferential loss of 
hydrogen could help the envelope remain bound. The 
planet evolution tracks presented do not include mass- 
loss. 

There do exist Neptune-size equilibrium configurations 
at T eq > 500 K for planets with H/He envelopes from 
core nucleated accretion. Our equilibrium planet struc- 
ture models in Section [2] explore and map out the (M env , 
M c , T eq , Ti n t) parameter space that yield radii within the 
range 2-6 i?© (Figure It is important to note, how- 
ever, that the equilibrium models provide "snap shots" 
of possible equilibrium configurations of planets under- 
going quasi-static evolution. The models do not address 
how a planet could reach a given state, nor the timescale 
for the planet to evolve out of the state. 

6.1.2. Outgassing of H2 

The second formation pathway to low-density 
Neptune-size planets we considered was outgassing of 
H2 from rocky planets. Outgassed low-mass planets 
(M p < 30 M©) without substantial H2O envelopes, how- 
ever, can only account for radii up to ~ 3 i?©. Even 
achieving 3 i?© with an outgassed envelope is a stretch, 
requiring (concurrently) a near-optimal initial planetesi- 
mal composition, full oxidization of accreted metals, and 
retention of most H2 released. Realistically, the major- 
ity of outgassed planets will be smaller than this radius 
limit, as we elaborate below. 

The metal and H2O content of the primordial mate- 
rial forming a planet set a strict limit on the amount 
of H2 that can be released via the outgassing reaction, 
2Fe + 3H 2 -> Fe 2 3 + 3H 2 . In this work, we have 
adopted a primordial chemical composition representa- 
tive of EH chondrites mixed with just enough addi- 
tional H2O ice to fully oxidize all the metals. Out of 
the Solar System chondrites, this composition should 
be near optimum for outgassing of H2 due to the 
high p roportion of unoxidized iron (in metal or sulfide 
form) (jElkins-Tanton fe Seagerll2008bl ). Typically plan- 
ets formi ng from a mixture of Solar System chondrite-like 
material (|Jaros cwich 1990, within which the proportion 
of metallic iron varies from 0.1 to 22 wt. %) would have 
a lower capacity to outgas H2. 

Even given a high initial amount of reduced metals, a 
planet's eventual outgassed envelope mass is contingent 
upon the fraction of metals that oxidizes. To bound the 
radii of outgassed planets, we considered the end-member 
case of complete oxidation of all Fe to Fe2 03. In this lim- 
iting case, the planet is core-les s with all its iron incor- 
porate d in the mantle as oxides dElkins-T anton fe Seagerl 
I2008a|) . Planets retaining a metallic core would degas less 
H2. Ultimately, the overall fraction of Fe that reacts with 
water is determined by the competition between the rate 
of oxidation and the rate of sinking of metallic Fe to form 
the planet iron core. For a more detailed discussion see 
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Fig. 9. — Energy limited mass loss rates for the planet models in Figure [4] Mass loss rates are estimated for the case where the planets 
orbit a star with similar properties to our sun (M* = 1 M Q , L* = 1 L Q , and L X uv = 3.4 x IO" 6 L ). A mass loss efficiency of e = 0.1 
is assumed. The line styles have the same meanings and correspond to the same model planets as in Figure [4] Each curve corresponds to 
a different value of M env /M p : 0.001 (thin solid), 0.01 (thin dashed), 0.05 (thin dot-dashed), 0.1 (thin dotted), 0.2 (thick solid), 0.3 (thick 



dashed), 0.4 (thick dot-dashed), and 0.5 (thick dotted). All data in this plot have L p /M p 
T eq = 1000 K. 



io- 



Wkg" 1 , and (a) T eq = 500 K or (b) 



lElkins-Tanton fc Seagerl (|2008al) . 

Finally, the mass-radius relations for outgassed planets 
in Section @] considered 100% retention of all outgassed 
H2. Atmospheric escape leads to less massive H2 en- 
velopes and smaller planets overall (Section [3]). Indeed, 
while the primary outgassed atmospheres surrounding 
Earth and Mars during their accre tion were likely H2- 
dominated (jSchaefer fc Fcglev 2013), both planets today 
harbor secondary atmospheres with higher mean molec- 
ular weights. 

How close can outgassed-planet radii plausibly get to 
the limiting outgassing M p — R p relation? Relaxing our 
assumptions of optimum outgassing conditions, we in- 
vestigate an intermediate, incomplete-oxidation case in 
which 50% of the accreted Fe is converted to FeO. This 
scenario leads to planets with 19.5% of their (initial) 
mass in an iron core, 0.6% by mass degassed in H2, 
and a mantle Mg# of 0.62 (Table EJ}. With no loss of 
H2 these planets could have radii up to 2.7 R® (again 
considering M p < 30 M®). Whereas, with atmospheric 
mass loss, planets that retain only 1 to 10% of the de- 
gassed H2 would have radii up to at most 2.4-2.5 i?® 
for T eq = 500 - 1000 K. Thus, radii up to - 2.5 R m are 
more realistically achieved by rocky planets with out- 
gassed H2 envelopes but no free water. Planets with a 
water layer between the rocky interior and H2 envelope 
could be slightly larger, but only if little or no H2 was 
mixed in with the H2O. 

6.2. Maximum Planet Radius at Specified Mass 

We have modeled the internal structure of low-mass, 
large-radius planets with hydrogen-dominated atmo- 
spheres. For planets with outgassed H envelopes, we 
derived a limiting low density M p — R p relation by lever- 
aging an upper bound on the amount of H2 that can be 
degassed from rocky planetesimals. The limiting low- 
density M p — R p relation is less clear cut for planets 
formed from core nucleated accretion, because the ini- 
tial reservoir of H/He accreted from the nebula need 
not be a constraining factor. Our detailed planet for- 



mation calculations provide discrete examples of planets 
at T eq = 1000 K with only a few Earth masses yet radii 
larger than Jupiter. 

The low-density limit for planets formed from core nu- 
cleated accretion depends on the heavy element core and 
envelope masses achievable at a given equilibrium tem- 
perature. The plausible combinations of {M env , M c , T eq ) 
in turn rely on the protoplanetary disk properties and the 
migration history of the planet. The heavy element core 
mass is determined by the isolation mass, given the solid 
surface density and the distance from the star where the 
planet forms. The isolation mass (and thus M c ) can have 
a wide range of values, from less than 0.1 M® to more 
than 20 M®. The initial mass of H/He accreted by the 
planet is determined by the availability of a gas supply 
from the disk as governed by disk lifetime relative to the 
time taken for the heavy element core to reach isolation 
mass. Disk lifetimes range over an order of magnitude 
- fro m 1 to 10 Myr, wit h a characteristic value of a few 
Myr (jHillenbrandi 120081 ) - leading to some freedom in 
the initial M env expected from core nucleated accretion. 
Mass loss over the planet's history would serve to de- 
crease M env over time. Finally, the current equilibrium 
temperature T eq depends on the migration history of the 
planet, and can, in principle, be anywhere from 100 K to 
2000 K. Thus, due to the large spread in observed disk 
properties, a wide range of (M env , M c , T eq ) from core nu- 
cleated accretion are plausible. We have shown detailed 
planet formation calculations for four reasonable choices 
of disk planetesimal densities and lifetimes. 

We have succeeded in placing a tighter constraint on 
the low-density M p — R p relation for outgassed planets 
than we have for planets from core nucleated accretion. 
This is due to the inherent limits on outgassed envelope 
masses; at very most, only a few percent of the mass of a 
planet can be outgassed in H2 . The end- member case of 
a planet that accreted from an optimum mixture of EH 
material and H2O ice, where all the water and iron re- 
acted, and where all released H2 was retained, sets an up- 
per bound on the transit radius possible at a given mass 
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for a rocky planet with out-gassed H2 atmosphere (Fig- 
ure [7j) . Typically, rocky planets with out-gassed H 2 at- 
mospheres would have mean densities above this limiting 
Alp — R p relation. It should be noted, that our limiting 
M p — R p relation applies to planets formed from material 
similar to Solar System chondrites. Planets formed from 
material with higher metallic Fe content would have the 
potential to outgas more H 2 . 

We have so far considered either core nucleated accre- 
tion or outgassing due to water-iron reactions as sepa- 
rate pathways for planets to acquire hydrogen rich en- 
velopes. Core nucleated accretion contributes near solar 
composition material to the envelope (Y ~ 0.25), while 
water-iron reactions contribute hydrogen but not helium 
(Y = 0). If both processes occur on the same exoplanet, 
an envelope with intermediate, sub-solar, non-zero He- 
lium content (0 < Y < 0.25) may result. 

The assumed chemical make-up of the planet envelope 
and heavy element core affect the planet M p — R p rela- 
tions for planets formed by core nucleated accretion and 
by outgassing. H/He envelopes in which He is depleted 
relative to solar will be more voluminous, for the same 
envelope masses, temperatures, and heavy element core 
properties. This is largely due to the influence of the 
mean molecular weight on the atmospheric scale height. 
For instance, decreasing Y=0.25 to Y=0.0 in the equi- 
librium planet models of Section [3.3[ increases the radial 
extent of the envelopes by ~ 15 - 20% for M p > 20 M®. 
For lower mass planets, the change in the gravitational 
acceleration between the top and bottom of the enve- 
lope can be substantial and Y can have a larger effect 
on the envelope thickness. Pure H envelopes can be up 
to twice as thick as the corresponding H/He envelope, 
near the low mass extreme of the M p — R p relations in 
Section r3.3l In our planet structure models, however, the 
effect of the envelope He abundance is largely offset by 
the higher density heavy element core composition in our 
outgassing models (EH chondrite cores) as compared to 
our core nucleated accretion models (ice-rock cores). 

We have mapped out the contribution of low-mass 
planets with hydrogen-dominated atmospheres to the 
limiting low-density M p (Rp) relation. Although we 
have not considered them in detail here, planets may 
also form with high molecular weight envelopes, for in- 
stance, after having accrete d large amounts of ices be- 
yond the snow line (e.g., iKuchnerl 120031 : iLeger et"aTI 
1200-11) . Higher molecular weight envelopes are more 
dense (with smaller atmospheric scale heights) than their 
hydrogen-dominated counterparts, but may be less af- 
fected by atmospheric escape. It is possible that planets 
with high molecular weight atmospheres could also con- 
tribute to the limiting low-density M p (R p ) relation for 
Neptune-size planets. 

6.3. Minimum Planet Mass at Specified Radius 

Our ideal goal was to determine a lower bound on the 
plausible planet mass given a planet radius in the range 
2-6 i?® and equilibrium temperature T > 500 K. We 
note that the relation defining the the maximum radius 
for a given planet mass does not necessarily translate 
into a relation for the minimum planet mass at a given 
radius. Indeed, at low masses, dR p /dM p < in the 
iso-composition M p — R„ relations for planets with gas 
envelopes (e.g., Figures |4] [5] and [7J) . Thus, in order to 



bracket the minimum planet mass of a transiting planet 
candidate, we must assess the survivability of low-mass 
planets for a range of interior compositions. 

Mass loss is a major limiting factor that constrains the 
minimum M p (Rp) for strongly irradiated (T > 500 K) 
Neptune-size planets harboring hydrogen dominated en- 
velopes (Section [5]) . This is true whether the planet ac- 
quired its envelope through core nucleated accretion or 
through outgassing. If the heavy element core mass is 
small (< 2 M®) and T eq is high (1000 K) then the planet 
will not be able to hold on to very much gas. With the en- 
ergy limited mass loss rates from Equation (|13j) . we may 
roughly assess the survivability of potential planet config- 
urations. By choosing a threshold envelope loss timescale 
^mo> we can derive a lower bound on the planet mass at 
a given radius based on the requirement t^ > t^ I0 . To 
illustrate this approach, we adopt t^ IQ = 1 Gyr and ex- 
plore what this implies for planets with low mean molec- 
ular weight envelopes from core nucleated accretion (Fig- 
ure [6]) and from outgassing (Figure [8]). 

We estimate, using Figure [6j the minimum masses 
of R p = 2-6 .R® planets with H/He envelopes formed 
by core nucleated accretion beyond the snow line. For 
R p = 6 i2® planets, the least massive planet models that 
satisfy t^ > 1 Gyr are 1.3 to 1.7 M® at T eq = 500 K, 
and 4.0 to 4.7 M® at T eq = 1000 K (for L p /M p between 
10~ 9 and 10~ n Wkg" 1 ). Analogously, at Rp = 4 i?®, 
the tfa > 1 Gyr survivability constraint requires that 
Mp > 1.1 to 1.4 M® at T eq = 500 K, and M p > 3.6 
to 4.3 M® at T eq = 1000 K. At Rp = 2 R m , almost 
all possible (M p , M env , L p , T eq ) configurations in Fig- 
ure [6] have sub-gigayear envelope-loss timescales, due to 
the small planet and envelope masses (M p < 5 M®, < 
Menv < 0.1% M p ). An ice/rock core surrounded by an 
H/Hc envelope from core nucleated accretion may not be 
a plausible interior composition scenario for 2 i?® planets 
at these equilibrium temperatures. Instead, other possi- 
bilities not considered here (e.g., high molecular weight 
envelopes, or envelope-less planets) may account for the 
minimum plausible planet mass at 2 _R®. 

We turn now to planets with outgassed hydrogen en- 
velopes but no surface water, and apply the envelope 
mass loss threshold to Figure |5J In addition to atmo- 
spheric escape, hydrogen-rich envelopes acquired by out- 
gassing are also constrained by the limited H 2 reservoir 
(magenta line in FigureE]). We find that, at T eq = 500 K, 
there exist potential planet configurations that satisfy 
*M — 1 Gyr with masses as low as 1 Af® for planet 
radii ranging from 2 to 3 i?®. Granted, these minimum- 
mass outgassing scenarios necessitate near maximal re- 
lease and retention of H 2 . In contrast, at T eq — 1000 K, 
all possible H2 envelopes leading to R p = 2 i?® have sub- 
gigayear envelope-loss timescales. For larger radii (2.5, 
and 3 i?®), planets with masses as low as 3.5 - 4 M® (de- 
pending on L p ) may pass the > 1 Gyr survivability 
criterion. 

We emphasize that minimum masses estimated follow- 
ing the approach above are contingent upon the chosen 
t^j threshold, the energy-limited mass loss parameter 
values assumed (here we took eLxuv / Lbol = 10~ 6 ), 
and the range of planet ages/intrinsic luminosities under 
consideration. High (lower) t^ I0 would leach to higher 
(lower) minimum M p (R p ). Although quantitatively very 
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assumption-dependent, minimum masses derived from 
may nonetheless yield important qualitative insights. 

6.4. Implications for Kepler Planet Candidates 

We conclude with a discussion of the implications of 
our results for the Neptune-size planet candidates dis- 
covered by Kepler. Candidates in the 2 — 6 i?® size-range 
account for a large fraction of the curren t candidates de- 
tected by Kepler (|Borucki et al.l l2011allbf ). This raises 
the question of why Neptune-size planet candidates are 
so common. One possible contributing factor revealed 
by this study is that not very much mass is needed in a 
hydrogen dominated envelope for a rocky heavy element 
core to reach radii within 2-6 i?® . 

Our main conclusion is that the Neptune-size planet 
candidates could have low mass (M p < 4 M®). This de- 
duction is supported by our calculations of the forma- 
tion, structure, and survival of planets with voluminous 
envelopes of light gasses. 

• Formation: We demonstrated that planets 3 to 
8 M® with substantial H/He envelopes can plau- 
sibly form by core nucleated accretion beyond the 
snow line and migrate to T eq ~ 500 K given rea- 
sonable disk surface densities and disk dissipation 
timescales. Migration to T eq ~ 1000 K with the 
envelope intact in timescales of a few Myr is more 
challenging. 

• Structure: We mapped the regions of (M p , M env , 
T eq , L p ) parameter space that yield radii between 
2 and 6 i?® for planets with H/He envelopes from 
core nucleated acretion and for planets with out- 
gassed H 2 envelopes (Figures[B]and[5J respectively). 
Since at most a few percent of a planet's mass 
can be degassed as H 2 , rocky super-Earths (M p < 



30 M®) with out gassed hydrogen atmospheres but 
without substantial H 2 typically will not account 
for Kepler planet candidates larger than ~ 3 i?®. 

• Survival: Envelope mass loss plays a major 
role governing the minimum plausible M p (R p ) for 
strongly irradiated (T > 500 K) Neptune-size plan- 
ets with hydrogen-dominated envelopes. At R p = 
2 f?®, H/He envelopes surrounding ice- rock cores 
would likely be lost in short order. At larger radii 
(2.5 to 6 i?®), planet configurations with envelope 
mass loss timescales longer than a Gyr (assum- 
ing eLxuv I Lbol = 10~ 6 ) exist down to masses 



1 M® at T ei 
T eq = 1000 K. 



500 K and down to ~ 4 M ffi at 



Neptune-size planets with masses M p < 4 M® could 
prove a challenge for radial velocity (RV) follow-up 
due to their low RV semi-amplitudes, but confirmation 
and mass measurements through transit timing varia- 
tions may be poss ible in some cases (e.g., Kepler-11, 
iLissauer et al.ll201~ll) . Figures [5] and [5] may be useful tools 
for assessing minimum masses for Kepler planet candi- 
dates. 
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